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Abstract 

We study the static, analytical solution of black holes in the warped DGP braneworld 
scenario. We show that the linearized field equations and matching conditions lead to 
solutions that are not compatible with Schwarzschild-(A)dS(4) solutions on the brane. 
This incompatibility is similar to vDVZ discontinuity in massive gravity theory. Fol- 
lowing the standard procedure to remove this discontinuity, which firstly was proposed 
by Vainshtein, we keep some appropriate nonlinear terms in the field equations. This 
strategy has its origin in the fact that the spatial extrinsic curvature of the brane plays 
a crucial role in the nonlinear nature of the solutions and also in recovering the well- 
measured predictions of General Relativity (GR) at small scales. Using this feature, 
we obtained an interesting black string solution in the bulk when it is compatible with 
4D GR solutions on the brane. 
PACS: 04.70.-s, 04.50.-h 
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1 Introduction 



A model universe with possible extra dimensions has attracted a lot of attention in recent 
years. Questions such as: are we living really in a four dimensional universe, or are there 
extra dimensions we cannot see, have not been addressed precisely yet. But, theoretically 
there are a lot of interesting outcomes by incorporating some extra dimensions in theories 
such as gravitation. In fact, extra dimensions are an integral part of fundamental theories of 
physics such as superstring/M-theories that need more than four spacetime dimensions. The 
idea to add an extra dimension to General Relativity's four dimensional spacetime, firstly was 
introduced by Kaluza and Klein to unify Electromagnetic and Gravity. Now p-branes and 
D-branes are well-known in the context of string theory where branes are solitonic solutions 
of 10-dimensional string theories, p-branes are extended objects of higher dimension than 
strings (1-branes) and the D-branes are kind of p-branes on which open strings can end [1]. 
In braneworld scenarios we assume that our (1 + 3)-dimensional spacetime is a domain wall 
embedded in a 5-dimensional spacetime called the bulk [1]. All matter fields and other gauge 
bosons live on the brane but gravitons, which carry the gravitational interaction, can travel 
into the extra dimension. There are at least two basic reasons for investigating braneworld 
models: firstly for solving hierarchy problem, that is, a large difference between Planck and 
Electroweak scales, secondly for proposing a unified theory of all fundamental interactions 
in the nature. 

During the last two decades, different kinds of braneworld scenarios have been pro- 
posed such as the Horava-Witten model [2], Arkani Hamed-Dimopoulos-Dvali (ADD) model 
[3], Randall-Sundrum (RS) models [4] and the Dvali-Gabadadze-Porrati (DGP) model [5]. 
Among these models, one of which has attracted most regards in last decade, is the RS one 
brane (RSII) model. In this model it is assumed that our universe is a 3-brane with posi- 
tive tension, the extra dimension is large and the bulk has a negative cosmological constant 
leading to a warped geometry. Another model that has attractive results from cosmological 
viewpoint, is the DGP model with one large extra dimension but the bulk is Minkowski. In 
this model, there is a new term in the total action that comes from quantum interaction 
between matter confined on the brane and the bulk gravitons which induces gravity on the 
brane. In this model, gravity leaks off the 4D brane into the bulk and becomes 5-dimensional 
at large scales, r 3> r c , but on small scales, gravity is effectively bounded to the brane and 4D 
dynamics is regained. The transition from 4D to 5D behavior is governed by the crossover 
scale r c . An interesting property of this model is that it contains a self-accelerating branch 
of the solutions which can explain late time acceleration due to a weakening of gravity at 
low energies [8]. 

We note that a basic requirement for any alternative theory of gravity, such as braneworld 
scenarios, is that they should reproduce General Relativity predictions in the appropriate 
limit to be phenomenologically viable. Two of the most important applications of GR are 
cosmology and black holes. Here, our concentration is on the issue of braneworld black holes, 
i.e. finding the bulk and the brane metric when a spherically symmetric energy-momentum 
distribution is localized on the brane. To obtain black hole solutions in a braneworld sce- 
nario, we note that generally braneworld solutions can be obtained by following two different 
approaches: In the first approach, dynamics and geometry of the whole bulk spacetime are 
primarily considered, then the dynamics on the brane is extracted by using the Darmois- 
Israel matching conditions. The second approach is to obtain effective 4-dimensional field 
equations on the brane firstly and then try to extent these solutions to the bulk. This ap- 
proach has been applied to braneworlds with induced gravity in Ref . [6] . The main difficulty 
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in the later approach is that it is not always possible to obtain a closed set of equations on 
the brane. So, we are not able to predict the behavior of the fields on the brane just with 
data on the brane. 



In 1999, Chamblin, Hawking and Reall found the Schwarzschild black string solution for 
RSII model as follows [10] 



In this solution, there is a line singularity along r = for all y, but the black string is unstable 
in response to large-scale perturbations because the 5D curvature is unbounded at the Cauchy 
horizon, as y — > oo. Another solution for static uncharged black hole was found by Dadhich 
et al. from the induced field equations on the brane [11]. In this work, a Reissner-Nordstrom 
type correction to the Schwarzschild solution was found and 5D gravitational effects, which 
are impressed by the bulk Weyl tensor, induce a tidal charge parameter Q. In this approach 
the bulk metric has not been found. Braneworld black holes have been studied extensively in 
recent years, some of these studies can be found in Refs. [10]-[17]. An extension of the DGP 
setup is the warped DGP scenario, which is a hybrid braneworld model containing both 
RSII and the DGP models as its limits [6, 7]. This model is constructed by considering the 
effect of an induced gravity term in the RSII model. Black hole solutions in a warped DGP 
braneworld have been studied recently in the case of a conformally flat bulk for spherically 
symmetric vacuum on the brane [18]. This study confirms the idea that an extra term in 
the effective vacuum field equations on the brane can play the role of a positive cosmological 
constant. 

Within this streamline, in this paper we present an analytical solution of static black hole 
in the warped DGP braneworld. Firstly, we consider pure (A)dS^) bulk field equations with 
a general spherically symmetric 5D metric to achieve independent field equations that should 
be held in the bulk. We take into account the 3-brane effect through appropriate matching 
conditions that act as boundary conditions on 3-brane universe in order to constraint general 
bulk solutions. In this way, we obtain a complicate system of equations which should be 
solved to reach the warped DGP braneworld, static, black hole solutions. Due to difficulties 
in solving the field equations analytically, we consider the weak field limit of the scenario 
to find those solutions of the field equations that describe the gravitational interactions in 
regions far enough from the gravitational source. We demand these solutions to support 
4D General Relativity solutions for r <C r c on the brane. However, the full linearized 
theory that we have adopted, does not lead to a correct 4D Schwarzschild- (A) dS( 4 ) solution 
on the brane. The reason for this incompatibility and its solution are summarized in which 
follows: The linear analysis of the DGP model shows that the tensor structure of the induced 
metric perturbations takes the five dimensional form even at short distances, i.e. in the limit 
r c — > oo, we have an incompatibility with Einstein gravity [5, 8, 9]. The situation is analogous 
to the case of models with massive gravitons in which we have a deviation from 4D GR 
(which is a massless theory), even in the massless limit. This discontinuity in graviton mass 
is known as van Dam-Veltman-Zakharov (vDVZ) discontinuity [19]. Vainshtein claimed that 
the vDVZ discontinuity is an artifact of the Pauli-Fierz Lagrangian, i.e. of the linearization 
of the true, covariant, non-linear equations of massive gravity, and the 4D GR is recovered 
by non-linear effect [20]. There have been many discussions about this issue, see for instance 
[21]. Another possibility that can disappear the mentioned discontinuity is introducing a 
cosmological constant [22]. In this way, in order to 4D GR be recovered, the cosmological 
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constant should be large enough as far as the graviton mass is not completely negligibly 
small. The paradox in DGP gravity seems to be that while it is clear that a perturbative, 
vDVZ-like discontinuity occurs in the potential somewhere, no such discontinuity appears 
in the cosmological solutions; at large Hubble scales, the theory on the brane appears safely 
like general relativity [8, 9]. Indeed, the cosmological solutions at large Hubble scales are 
extremely nonlinear, and that perhaps, just as Vainshtein suggested for massive gravity, 
nonlinear effects become important in resolving the DGP version of the vDVZ discontinuity. 
In the context of DGP model, there are various works that indicate 4D GR is recovered at 
short distances [23, 24]. It was shown that the brane bending becomes non-linear at scale 
shorter than r* = (r^) 1 / 3 , thus the linear analysis breaks down there. Also the leading 
order correction to 4D GR has been obtained in this regime [24]. In Ref. [25], an alternative 
formalism which can handle general perturbations in weak gravity regime was proposed. 
They generalized the DGP model by taking into account the bulk cosmological constant 
and the brane tension balanced with it, so the background metric is AdS( 5 ) in the bulk 
and Minkowski on the brane. By perturbing the background metric and following the prior 
approach that had been proposed for RS scenario [26], they confirmed the recovery of 4D 
GR at short distances and re-derived the leading order correction to it. They have done this 
work by taking into account the non-linear brane bending for weak gravity at small scale 
r < r*. 

Here, in order to recover 4D Schwarzschild-(A)dS(4) solution on the brane for r < r* or in 
the limit r c — > oo, with regard to previous works, we take into account some nonlinear terms 
in the field equations. In fact, these nonlinear terms become so large in that regions that 
one cannot neglect them relative to other nonlinear terms. We show that nonlinear effects 
become important in resolving the "Warped DGP" version of the vDVZ discontinuity too. 
By considering these nonlinear terms, we have found interesting solutions for static Black 
holes on the warped DGP brane. 

2 The field equations 

We consider a 3-brane £ embedded in a 5-dimensional bulk M.. The total action for the 
system is 

S = A / <FXV=9 ( 5) R - 2A 5 ) + / d 5 X^ C bulk + S hrane , (1) 
where S^. ane is the 3-brane action defined as 

Sbrane = 77^ I d^x^q ( ^ - 2A 4 ) + f d^x^q C brane + f d^x^q ^ . (2) 

cjab is the 5D bulk metric with corresponding Ricci tensor given by ^Rab- The brane has 
induced metric q^ u with corresponding Ricci tensor ^R^ u . The boundary Gibbons-Hawking 
term is implied to yield the correct Einstein equations in the bulk. A5 and A4 are bulk 
and brane cosmological constants respectively. % can be interpreted as the brane tension 

of the standard Dirac-Nambu-Goto action and can include quantum contributions to the 
four-dimensional cosmological constant [12]. Cbulk and Cb rane are the bulk and brane matter 
Lagrangian respectively. We note also that the Planck masses M 5 and M4 are defined as 
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k\ = 87tG( 5 ) = M 5 3 and k\ = 8nG^) = M 4 2 , where M 5 and M 4 have dimensions of 
(length) -1 . 

We chose a coordinate y for extra dimension so that our 3-brane in the Gaussian Normal 
Coordinates is localized at y — 0. By variation of the action with respect to the bulk metric, 
5D field equations would be obtained as follows 

^G AB = -A 5 g AB + 4 ^T AB + 5(y)4 ^T AB , (3) 

where 

{loc) T AB = g\g B ( - 1) ^ g ( 4) G,u + A 4 ^ - 4 ^ ane % u ) , (4) 
is the localized energy-momentum tensor on the brane and 



^T AB = -2^f + g AB C bulk . (5) 



^G AB and ^G^ u denote the Einstein tensors constructed from the bulk and the brane 
metrics respectively. The tensor q^ u is the induced metric on the brane E. The field equations 
in the bulk take the following form 

^G AB =^R AB - l - MRgju = -A 5 g AB + 4 ^T AB , (6) 

or 

(S) R B _ 1 <5) R5 B = _ A5§ B + K 2 H) r fl (7) 

Moreover, the following modified (due to the presence of induced gravity on the brane) 
Israel-Darmois junction conditions would be obtained from the distributional character of 
the brane content 

2 

[K, v ] - q^[K] = -4 ^T, v = % ( (4) GV + A 4 g J - 4 ^ brane % u , (8) 

K 4 

where K^ u = \d y {g^ v ) is the extrinsic curvature of the brane and brackets denote jump 
across the brane (y — 0). 

Now we consider a general, 5-dimensional, static metric with spherical symmetry on the 
brane as follows 

ds h 2 = -e v ^ y) dt 2 + e x{r ' y) dr 2 + r (d6 2 + sin 2 6d(j) 2 ) + dy 2 . (9) 

Also we assume the matter content of the bulk is just a cosmological constant A 5 (i.e., 
(buik)rp AB _ g^j anc j foe ma ^ er content of the 3-brane universe is considered to be a cosmo- 
logical constant A 4 plus a localized spherically symmetric untilted perfect fluid (e.g. a star) 
with 

(brane)^ = (p + + ^ ^ (1Q) 

where stands for the 4-velocity of the fluid and p = p = for r > R. Nevertheless, since 
we want to obtain static black hole solutions outside the star ( that is, for r > R), in these 
regions the brane has just a cosmological constant A 4 . 



5 



Hereafter, we consider a Z2-symmetry on reflection across the brane, thus at r > R the 
Israel matching conditions become 

KZ-K5Z = r e ( A) Gr+AtS»), (11) 

where r c = ^ = is the DGP crossover distance, and by definition = K^(y = + ) = 

—K^(y = CT). Now by putting the metric (9) into the bulk field equations (7) we obtain 
the (tt), (rr), (00) and (ry) components of the bulk field equations respectively as follows 

4r/i rr + 3rfi 2 r + 12^ r - 4A r - 2rX r ^ r + -(l - e A ~ M ) + 
2r\ yy e x + Arfi yy e x + r\ 2 y e x + 3rfi 2 y e x + 2r\ y fi y e x = 4rA 5 e A , (12) 

4 . . 

r/4 + 2rv r /j, r + Au r + Afi r + - (1 - e A ~^J + 2rv yy e x + 
Arfi yy e x + ru 2 e x + 3r^e A + 2rv y n y e x = 4rA 5 e A , (13) 

2rv rr + 2rfi rr + rv 2 + r/i 2 . — ru r X r + rv r ^i r — rA r /i r + 2v r — 2\ r + 4/i r + 
2r\v yy + X yy + ix yy \e x + r\\ y v y + \ yi x y + v yi i y \e x + r\\\ + + ix 2 y \e x = 4rA 5 e A , (14) 

4rii rr + 2rv rr + 3r^i 2 + rv 2 + 2rz/ r /i r — 2r\ r fi r — ru r \ r + 4z/ r — 4A r + 12/i r + 
4 / \ 

- (l - e A ^J + rfx 2 y e x + r[\ y v y + 2\ y ^i y + 4i/„^]e A = 4rA 5 e A , (15) 

2z/ r?/ + Afx ry + - v r \ y + 2^ r ^ - 2A ?/ /i r + - (jiy — A y J = . (16) 

By using the Israel matching conditions (11), we obtain the following conditions on the 
brane 

1 6 A 

— -(2^ y +A y )| J/=0 + = r c {A 4 - — {A-Ae x -^+4r 2 fi rr +3r 2 fi 2 -4r\ r +12rfi r -2r 2 X r ^ r )} , (17) 

1 e A 

- - (2^ + | y=0 + = r c {A 4 - ^ (4 - 4e A ^ + r 2 ^ + 4r// r + 4rz/ r + 2r V/v) } , (18) 

1 e A 

~2 { Uy + ^ + Vv) \y= 0+ = r c{ A 4 - — {2rv rr + 2rfi rr + r^i 2 . + rv 2 . 

—rv r \ r + 2v r — 2A r + 4/i r + ru r jj, r — r\ r jj, r ^j j . (19) 

The subscripts y and r in these relations represent partial differentiation with respect to 
y and r respectively. Note that these equations are held on the brane outside our spherical 
object. 
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3 The weak field solutions 



The nonlinear, second order partial differential equations (12) -(16) should be solved and 
the corresponding solutions should satisfy junction conditions (17), (18) and (19) across the 
Z2-symmetric brane. This is not an easy task at all! To find some analytical solutions, 
we consider just the weak field regime (i.e., far enough from the source localized on the 
brane). In this respect, we adopt the assumption that |A| and \fj,\ are small quantities 
compared to unity; that is, |A|, \fi\ <C 1. By adopting this assumption, we linearize 
our field equations and also brane boundary conditions. By keeping only the leading order 
terms, the bulk equations (12) - (16) reduce to the following equations 

c 4 
2r/i rr + 4^ r - 2A r - v r + r(\ yy + 2/j, yy ) = -rA 5 , (20) 

T O 

c 4 

Ur ~ ~2 + r ^yy + = 3 rAs ' ( 21 ) 

r(v„ + fi rr ) + v r + 2/i r - A r + r(uyy + \ yy + fi yy ) = 2rA 5 , (22) 



-(A-/i) = v r + 2pi r + — - -rA 5 , (23) 



2 /x c 2 

-(A - /i) = i/ r + 2/v + — - -i 

r r l 3 

where c is an integration constant to be fixed later. In the same procedure applied to 
equations (17), (18) and (19), we achieve the following linearized junction conditions 

1 

~2 



(X y + 2(j,y)\ y=0 + = r c |A 4 - — (ji - A + 3r/i r + r 2 /i rr - rA r )|, (24) 
-^(^ + 2^)1^=0+ =r c {A 4 --^( / [x-A + r / u r + rz/ r )}, (25) 

_ 2^f + ^ + A*2/)l2/=0+ = ^c{A 4 - 7^( rz/ rr + J~/i rr + 2^ r + l/ r - A r )}. (26) 

By solving the bulk field equations (20) - (23), we find \x and A in terms of v as follows 



\t{r,y) = - — J J iy r dydy-^u+^y 2 + ^A 5 y 2 +^F 1 (r) + H(r^j\y\+F 2 (r) + G(r) , (27) 

A(r, y) = -2 J J v„dydy ~-fj "rdydy - 2v + + ^A 5 r 2 + r\y\F[(r) + 

W^r) + rF 2 '(r) + F 2 (r) + ^A 5 y 2 - ^y 2 + A(r)\y\ + B(r) , (28) 

where A(r), B(r), Fi(r), F 2 (r), G(r) and H(r) are arbitrary functions of r, and a prime 
denotes derivative with respect to r. Also u(r, y) should satisfy the following differential 
equation 

-— -I - — 
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Therefore, for a given v(r,y) that satisfies this equation, //(r, y) and A(r, y) can be ob- 
tained via (27) and (28) respectively. Equation (29) has the following solution for i/(r, y) 

z/(r, y) = — + br 2 + ey 2 + d\y\ + — \y\ + qr 2 \y\ . (30) 
Now by using equations (27) and (28), we calculate \(r,y) and fi(r,y) to find 
w x / f 5 \ , ,o A a + c) b 4 , \ 9 

(2 9 r 2 + -/ + 2d- rf; - J\ - P(r)) + - 2tr 2 + f A 5 r 2 + rfj + F 2 + Q(r) , (31) 

V r ' ).r A 



2r 3 
/ \ ( f Q\< is /( a + c ) 1 * & 



y 2 - 



+ f r 2 + ^ " Fi " Af (r)) H - £ - V + F 2 + N(r) , (32) 

where a, b, d, e, / and q are constants, whereas -P(r), Q(r), M(r) and 7V(r) are arbitrary 
functions of r. Equations (20) - (26) can help us to determine the arbitrary parameters and 
functions appeared in u(r,y), X(r,y) and fi(r,y). By this way, we found that b, d, f and 
q should be zero and also A 5 = — A 4 and e = ^ A 5 in order for equations (20) - (26) to be 
satisfied. In this situation the solutions simplify to the following equations 

u(r,y) = ^ + ^A 5 y 2 , (33) 

A(r, y) = -(^> 2 + r c (^A 5 + ^±£) \y\ + ^ - ^A 5 r 2 + rW + W(r), (34) 

. , /a + c\ 2 /4 a + c\, a . orN 

Mr, v) = (^Jl/ 2 + r c (g A 5 - -^r) |y| + W(r) - - , (35) 

where W(r) = A r (r) + F 2 (r). The corresponding quantities on the brane (y = 0) take the 
following forms 

u(r,y = 0) = -, (36) 
r 

A(r, y = 0) = ^ + ^A 4 r 2 + rW + W(r), (37) 

rir,y = 0) = W(r)-£. (38) 
Since we expect /z(r, y = 0) = 0, we set W(r) = f-. So we find 

X(r,y = 0) = C ^+ 1 -A 4 r 2 . (39) 

c is an integration constant which we set to be zero in order to find more familiar results 
as follows 
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is(r,y) 



~ + 7tA 5 ?/ 2 

r 3 



(40) 



\(r,y) 



a l.o a 9 /4, a \ . 



(41) 



(42) 



Note that in these relations we have used A 5 = — A 4 interchangeably. To find Schwarzschild 
solution in appropriate limit, we set a=—2m , where m is the Schwarzschild geometric mass. 
It is clear that in our linearized warped DGP setup, A(r, y — 0) ^ y = 0) on the 

3-brane. We expect that in the limit of r c — > oo the solutions reduce to the 4D general 
relativistic solutions on the brane, that is, to a Schwarzschild- (A) dS( 4 ) metric. However, 
this has not occurred here. We note that the reason behind having a non-Schwarzschild- 
(A)dS( 4 ) solution can be explained through the fact that the curvatures of the bulk and brane 
spacetimes should not exactly cancel each other. In the next section we focus on this issue. 

4 Recovering the Schwarzschild solution on the brane 

To recover the Schwarzschild- (A) dS (4) solution on the brane, we note that all metric compo- 
nents should be much smaller than unity as a result of the weak field prescription. But this 
is violated by the linearized solutions we obtained in the last section. Nevertheless, if we 
look at the 5D bulk solutions Eqs. (40), (41) and (42), we will see that in 5D, A and \i are 
large for large r c , ( see for instance [24], [23] and [9]). Although we are interested only in 4D 
solutions, the linearized theory is not applicable at large r c . These features show that the 
completely linearized theory does not result in correct 4D Schwarzschild-(A)dS( 4 ) solutions 
(which are given by /i(r) = 0, A(r) = —v(r) and u(r) = — ^ — lA 4 r 2 ). In fact, in the 
limit r c — >oo, but yet in the weak field regime, we cannot neglect all nonlinear terms in some 
subspaces. For instance, if r\\a\ ^> r 3 , quantities such as A^ and fi y are large enough that 

we have to save their squares in the field equations. Therefore, for \a\ <C r <^ r* = (r^|a|)3 
the nonlinear field equations that should be solved are as follows* 



4r/i rr + 8fi r - 4A r - 2v r + 2U{r) + 2r(X yy + 2fx yy ) + r\ 2 y + 3r^ + 2r\ y fi y = 4rA 5 , (44) 



*We note that a static source placed on a brane creates a nonzero scalar curvature around it. For a source 
of the size < , this curvature extends to a distance ~ r* . More intuitively, a static source distorts a brane 
medium around it creating a potential well, and the distortion extends to a distance r <~ r* . This curvature 
suppresses nonlinear interactions that otherwise would become strong at the scale below r* , see for instance 




(43) 



2v r + 2U(r) + 2rv yy + Ar/i yy + 3r/£ = 4rA, 



(45) 



[14]. 



9 



2rv rr + 2rfi rr + 2v r - 2A r + 4/i r + 2r(z/ ra + A ra + /%,) + r\ y n y + r(\ 2 + fx y ) = 4rA 5 , (46) 

Arfi rr + 2rv„. + 2z/ r - 4A r + 8fi r + 2U (r) + 2r A^ + rji 2 y = 4rA 5 , (47) 

where U (r) is an arbitrary function of r . Junction conditions in this situation are the same 
as the linearized case given by Eqs. (24), (25) and (26). By solving this set of equations, we 
find the following class of solutions 

u(r, y) = + ^A 5 r 2 + ^r c A 5 \y\ + ^ 2 c A 2 5 y 2 , (48) 
a 1 4 

A(r,y) = -A 5 r 2 + -r c A 5 |y| , (49) 

v y y 

4 

A*( r >f) = g r cA 5 |j/| , (50) 

where now A4 = — IA5 with A4 = Note that in this framework we obtained a 

geometric interpretation of the brane cosmological constant in terms of the DGP crossover 
scale. Moreover, as we have mentioned at the end of the previous section, the curvatures of 
the bulk and brane spacetimes do not cancel each other exactly which is the reason behind 
having a non-Schwarzschild-(A)dS( 4 ) solution in the previous section. On the brane with 
y = 0, we find 

K^ = 0) = --iA 4 r 2 , (51) 
r 6 

X(r,y = 0) = -- + ]-A 4 r 2 , (52) 
r 6 

Mr,?/ = 0) = 0, (53) 

which is compatible with 4D Schwarzschild-AdS( 4 ) solution at the leading order if we set 
a = — 2m 

ds 2 = -(l - — - -A 4 r 2 ) dt 2 + -. — r + r 2 dVL 2 . (54) 

V r 3 ) (l-lf-lA^ 

Therefore, the Einstein gravity is recovered and the vDVZ problem is resolved by this 
strategy. Indeed, the spatial extrinsic curvatures of the brane, i.e. W y =o and ^ y \ y =Q , play a 
crucial role in the nonlinear nature of the solution and recovering the predictions of general 
relativity. In this respect, the nonlinear behavior arises from purely spatial geometric factors 
(see also Ref. [24]). We note also that there is a massless scalar mode, which arises from 
decomposition of five degrees of freedom of the bulk graviton in the massless limit (i.e. 
r c — t>oo). This extra scalar mode persists as an extra degree of freedom in all regimes of the 
theory. The extrinsic curvature suppresses this extra scalar field inside the region r<Cr t , 
whereas outside this region, the brane bending mode is free to propagate (see Ref. [9] for 
more details). 

We note that Eqs. (48), (49) and (50) can be rewritten by using A 4 = — |^ as follows 
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, , 2m lr 2 1 . 3 y 2 

^y) = -— + -^ + -\y\ + -^ c , (55) 

2m 1 r 2 1 . . . 

X ^y) = — -- A -2 + -\y\i (56) 

^,y) = — Is/I - (57) 

In the limit r c — >■ oo we find the Schwarzschild solution on the brane. This means that 
the effect of the brane cosmological constant is negligible in this limit. Thus, it is essentially 
impossible to observe a Schwarzschild- (A) dS (4) interaction in our observations, and the New- 
tonian gravity is dominant at least in our solar system scale. In is important to note that the 
above solution is actually a black string solution since we find a singularity at r = for any 
value of y. This can be shown through singularity of the 5-dimensional Kretschmann scalar. 
The expression of the 5-dimensional Kretschmann scalar is to long to be presented here, but 
our calculation by using the Maple package, has shown that the 5-dimensional Kretschmann 
scalar is singular at r = for all values of y. In fact, in the r — > limit, it reduces to the 
following simpler form 

g2= 32rg y +16,y + l^l 
r^o (r c + y) 4 r 4 

which obviously diverges at r = 0. This shows that singularity at r = is an intrinsic 
singularity, and this is true for all values of y, which confirm that the solutions are actually 
black string. Moreover, the 4-dimensional Kretschmann scalar of the metric (54) expressed 
in terms of r c is given by 

K^('^ s (^ = 3 - ^ + f m2 \ (59) 

which shows that singularity at r = on the brane is an intrinsic singularity. In the limit 
of r c — >oo, we have just one horizon which is the Schwarzschild horizon at r = 2m. But for 
m = there is no de Sitter horizon on the brane. The stability of black string solutions 
should be considered in the light of the results of Gregory-Laflamme pioneering work [27]. 

Another issue which requires especial attention here is the fact that a de Sitter bulk 
invalidates most of the attractive points of the Randall-Sundrum model such as the local- 
ization of the zero mode of the graviton near the brane. The question then arises: how is 
the situation in the warped DGP setup? In DGP-like models that have a 4D Ricci scalar in 
the bulk action, gravity becomes 5-dimensional at r ^> r c where 5D Einstein-Hilbert action 
is dominant. For r ^ r c gravity is 4-dimensional but not the 4D GR one. This property of 
these models is definitely different from RSII model, in which at low energies, Hi <C 1, 4D 
gravity is recovered to a good approximation. In RSII model what prevents gravity from 
leaking into the extra dimension at low energies is the negative bulk cosmological constant. 
Moreover, in RSII model the 4D gravitational constant strongly relies on the presence of 
the brane tension, and the brane tension should be positive to effective 4D gravitational 
constant have correct sign. This positive tension and the bulk negative cosmological con- 
stant are balanced with each other to have a zero 4D cosmological constant on the brane. 
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In our warped DGP setup, we found A4 < and A5 > 0, which is completely different from 
the RSII case. Here we don't need a negative A5 and positive A4 because induced gravity 
term gives us the desired results. Solutions that we have found in Sec 4, Eqs. (55)-(57), are 
valid only at r <C r* and for r ^> r c appropriate solutions should be obtained. Nevertheless, 
since a warped DGP scenario is a hybrid braneworld model which contains both the pure 
DGP and RSII models in appropriate limits [6], it is expected that some of the mentioned 
shortcomings are still present in this case too. This needs further justification that we are 
going to study separately. 

As the final remark, we note that Birkhoffs theorem is absent in theories of modified 
gravity such as the DGP braneworld scenario [28]. By calculation of the gravitational force 
on a test particle due to a spherical mass shell in the DGP setup, one can show that unlike 
in GR, the force depends on the mass distribution. In particular, the gravitational force 
within a spherical mass shell depends on the geometric structure of the bulk. In our setup, 
having Schwarzschild-(A)dS( 4 ) solution means that Birkhoff's theorem holds good and this 
is due to geometric structure of the bulk manifold in this case which differs from pure DGP 
case. 

5 Summary 

In this paper we obtained a class of static black hole solutions in the warped DGP braneworld. 
Firstly we solved the general bulk field equations in the weak field limit (by linearizing the 
field equations and matching conditions). The solutions on the brane then were found by 
using the Israel matching conditions. However, these brane solutions were not compatible 
with Schwarzschild-(A)dS( 4 ) solutions on the brane. We adopted a strategy to solve this 
problem based on keeping appropriate nonlinear terms in the field equations. This strategy 
has its origin in the fact that the spatial extrinsic curvature of the brane plays a crucial role 
in the nonlinear nature of the solutions, and also in recovering the predictions of General 
Relativity. In fact, the nonlinear behavior arises from purely spatial geometric factors. Using 
this feature, we obtained some interesting black sting solutions compatible with well-known 
black hole solutions on the brane. 
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